Effective action, magnetic excitations and quantum fluctuations in lightly doped 

single layer cuprates 
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We consider the extended 2D t — t' — t" — J model at zero temperature. Parameters of the model 
corresponds to doping by holes. Using the low doping effective action we demonstrate that the 
system can f) preserve the long range coUinear antiferromagnetic order, 2) lead to a spin spiral 
state (static or dynamic) , 3) lead to the phase separation instability. We show that at parameters of 
the effective action corresponding to the single layer cuprate La2-i:Sra;Cu04 the spin spiral ground 
state is realized. We derive properties of magnetic excitations and calculate quantum fluctuations. 
Quantum fluctuations destroy the static spin spiral at the critical doping Xc ~ O.f f. This is the 
point of the quantum phase transition to the spin-liquid state (dynamic spin spiral). The state is 
still double degenerate with respect to the direction of the dynamic spiral, so this is a "directional 
nematic". The superconducting pairing exists throughout the phase diagram and is not sensitive 
to the quantum phase transition. We also compare the calculated neutron scattering spectra with 
experimental data. 

PACS numbers: 74.72.Dn, 75.10.Jm, 75.50.Ee 



I. INTRODUCTION 

The phase diagram of the prototypical cuprate su- 
perconductor La2-xSr2;Cu04 (LSCO) shows that the 
magnetic state changes tremendously with Sr doping. 
The three-dimensional antiferromagnetic (AF) Neel or- 
der identified^ below 325 K in the parent compound dis- 
appears at doping x ~ 0.02 and gives way to the so- 
called spin-glass phase which extends up to x ~ 0.055. 
In both, the Neel and the spin-glass phase, the system 
essentially behaves as an Anderson insulator and exhibits 
only hopping conductivity. Superconductivity then sets 
in for doping x > 0.055, see Refj^. One of the most 
intriguing properties of LSCO is the static incommen- 
surate magnetic ordering observed at low temperature 
in elastic neutron scattering experiments. This ordering 
manifests itself as a scattering peak shifted with respect 
to the antiferromagnetic position. Very importantly, the 
incommensurate ordering is a generic feature of LSCO. 
According to experiments in the Neel phase, the incom- 
mensurability is almost doping independent and directed 
along the orthorhombic b axis^. In the spin-glass phase, 
the shift is also directed along the b axis, but scales lin- 
early with doping"*'^'^. Finally, in the underdoped su- 
perconducting region (0.055 ^ x < 0.12), the shift still 
scales linearly with doping, but it is directed along the 
crystal axes of the tetragonal lattice^. Very recent stud- 
ies reveal also the evolution of inelastic neutron spectra 
with dopingS.. 

Near x = 0.12 certain La-based materials develop a 
strongly enhanced static incommensurate magnetic order 
accompanied by small lattice deformation at the second 
order harmonic a^dOd^ , see also Ref for a review. 

Incommensurate features have also been observed 
in inelastic neutron scattering from YBa2Cu306+j^ 
(YBCO)i^iiiii^'i^Jiii^J^. In underdoped YBCO there 



is a rather large uncertainty in the determining of the 
doping level. However, it seems that the incommensu- 
rability in YBCO is about 30-40% smaller than that in 
LSCO comparing the same doping level. In a very recent 
work^° the electronic liquid crystal state in underdoped 
YBCO has been reported. The state has no static spins, 
but nevertheless, it demonstrates a degeneracy with re- 
spect to the direction of the dynamic spin structure. In 
addition, there are indications that the electronic liquid 
crystal state observed in^" is very close to a quantum 
phase transition to a state with static spins. 

The 2D t — J model was suggested two decades ago 
to describe the essential low-energy physics of high- Tc 
cuprates2iii2^i2^. In its extended version, this model in- 
cludes additional hopping matrix elements t' and t" to 
2nd and 3rd-nearest Cu neighbors. The Hamiltonian of 
the t — t' — t" — J model on the square Cu lattice has the 
form: 

H = —t — t c|^CjV — t cl^Cj"a 

{ij')<y 

+ J ^ (s.S, ~ jN^N,) . (1) 

Here, cj^. is the creation operator for an electron with spin 
a (cr =tji) at site i of the square lattice, (ij) indicates 
1st-, {ij') 2nd-, and {ij") 3rd-nearest neighbor sites. The 
spin operator is = ^cl^aapCip, and N.^ = J2a4a^io- 
with {Ni) = 1 — x being the number density operator. In 
addition to the Hamiltonian ([T|) there is the constraint of 
no double occupancy, which accounts for strong electron 
correlations. The values of the parameters of the Hamil- 
tonian ([1]) for LSCO are known from neutron scattering^, 
Raman spectroscopy^i and ab-initio calculations^^. The 
values are: J « 140 meV, t « 450 meV, t' « -70meV 
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, and t" « 35meV. Hereafter we set J = 1, hence we 
measure energies in units of J. 

The idea of spin spirals in the t — J model at finite 
doping was first suggested in Refi^. The idea had ini- 
tially attracted a lot of attention, see e. g. Refsi^s^i^. 
However, it has been soon realized that there was a 
fundamental unresolved theoretical problem of stabil- 
ity of the spiraP^. Together with lack of experimen- 
tal confirmations this was a very discouraging devel- 
opment. The observation of static and quasistatic in- 
commensurate peaks in neutron scattering caused a re- 
newal of theoretical interest in the idea of spin spirals in 
cuprates^Si^i^i^^iMi^^i^i^i^. It has been realized that 
in LSCO the charge disorder related to a random distri- 
bution of Sr ions plays a crucial role and in the insulating 
state, X < 0.055, the disorder qualitatively influences the 
problem of stability of the spiral. The point is that in the 
insulating state the mobile holes are not really mobile, 
they are trapped in shallow hydrogen-like bound states 
near Sr ions. The trapping leads to the diagonal spin 
spiral^'^'^'^. Percolation of the bound states gives 
way to superconductivity and in the percolated state the 
spin spiral must be directed along crystal axes of the 
tetragonal lattice^^,. So the percolation concentration is 
Xper = 0.055. The rotation of the direction of the spin 
spiral is dictated by the Pauli exclusion principle. The 
disorder at x > 0.055 is still pretty strong. However, un- 
like in the insulating phase, the disorder does not play a 
qualitative role and therefore in the first approximation 
one can disregard it. Thus, we arrive at the case of small 
uniform doping. This is the problem we address in the 
present work. 

As we already mentioned, the case of an uniform spin 
spiral (no external disorder) in a doped quantum antifer- 
romagnet has an inherent theoretical problem. If consid- 
ered in the semiclassical approximation, the out-of-plane 
magnon is marginal and in the end this implies an insta- 
bility of the spin spiral^^. An attempt to fix the prob- 
lem by account of quantum fiuctuations within the 1/S 
spin-wave theory was done in Ref^i. We understand now 
that, while being qualitatively correct, the work^ did not 
account for all relevant quantum fiuctuations. The effec- 
tive action method is much more powerful then the 1/S 
expansion because the method accounts for all symme- 
tries exactly and generates a regular expansion in powers 
of doping X, this is the true chiral perturbation theory. 
This is why in the present work we employ the effective 
action method. 

The structure of the paper is the following. In section 
in] we discuss the effective low-energy action of the modi- 
fied t-J model. Section UTTl addresses the issue of stability 
of the Neel state under doping. The spiral ground state 
in the mean-field approximation is considered in section 
ITVl The in-plane magnons are discussed in section IVl 
and out-of-plane magnons in section |VT1 Section |VlTl ad- 
dresses the quantum fluctuations and the quantum phase 
transition to the directional nematic. Finally discussion 
and comparison with experiments is presented in the sec- 



tion IVnTl 



II. EFFECTIVE LOW-ENERGY ACTION OF 2D 
t~t' ~t" - J MODEL AT SMALL DOPING 

At zero doping (no holes), the t-J model is equivalent 
to the Heisenberg model and describes the Mott insula- 
tor La2Cu04. The removal of a single electron from this 
Mott insulator, or in other words the injection of a hole, 
allows the charge carrier to propagate. Single-hole prop- 
erties of the t-J model are well understood, see Ref42i for 
a review. A calculation of the hole dispersion at values of 
parameters t, t' , and t" corresponding to the single layer 
cuprate LSCO has been performed in Ref4i using the 
Self Consistent Born Approximation (SCBA), see also 
Ref.'^^. According to this calculation the dispersion of 
the hole dressed by magnetic quantum fluctuations has 
minima at the nodal points qo — (±7r/2, ±7r/2), and it is 
practically isotropic in the vicinity of each point, 

^{P)^l(3p\ (2) 

where p = q — qo- We set the lattice spacing to unity, 
3.81 A 1. The SCBA approximation gives /3 « 2.2 « 
300 meV . The effective mass corresponding to this value 
is approximately twice the electron mass and this agrees 
with recent measurement of Shubnikov - de Haas oscil- 
lations^. In the present work we use /? as a fitting pa- 
rameter. We will see that to fit inelastic neutron data at 
a; = 0.1 we need 

(3^2.7. (3) 

This agrees well with the value obtained within the 
SCBA. The quasi-particle residue Z at the minimum of 
the dispersion is Z w 0.38'^^. In the full-pocket descrip- 
tion, where two half-pockets are shifted by the AF vector 
Qaf = (ti", tt), the two minima are located at Sa — -f) 
and Sb = (§' ~f ) • The system is thus somewhat similar 
to a two- valley semiconductor. 

The relevant energy scale for small uniform doping at 
zero temperature is of the order of xJ <^ J, relevant mo- 
menta are also small, q <^ I. Hence, one can simplify the 
Hamiltonian of the t-J model by integrating out all high- 
energy fiuctuations. This procedure leads to the effective 
Lagrangian or effective action. The effective Lagrangian 
has been first discussed quite some time ago ^^i"^^'"^^ , see 
also a recent work**-. That discussion resulted in the kine- 
matic structure of the effective Lagrangian valid in the 
static limit^^. This limit is sufficient only for the mean- 
field approximation. The time-dependent terms that are 
necessary for excitations and quantum fluctuations have 
been derived only recently^. The effective Lagrangian 
can be written in terms of the bosonic n-field that de- 
scribes the staggered component of the copper spins and 
in terms of fermionic holons ^. We use the term "holon" 
instead of "hole" because spin and charge are to some 
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extent separated, see discussion below. The holon has a 
pseudospin that originates from two sublattices, so the 
fermionic field V' is a spinor acting on pseudospin. For 
the hole-doped case, the effective Lagrangian reads 



(4) 



- i/'lea(7')^„ + %/2.g(V'laV'a)-["X (e„-V)n]| . 

The first two terms in the Lagrangian represent the usual 
nonlinear a model (NLSM), the field n is the subject 
of the constraint = f. The magnetic susceptibility 
and the spin stiffness are xi^ ~ 0.53/8 « 0.066 and 
Ps K, O.fS^. The rest of the Lagrangian in Eq. (jl]) rep- 
resents the fermionic holon field and its interaction with 
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the n-field. The coupling constant is 
The index a = a, 6 (flavor) indicates the location of 
the holon in momentum space (either in pocket Sa or 
Sb). The kinematic structure of the coupling term was 
first derived in Ref.— . The operator ct is a pseudospin 
that originates from the existence of two sublattices and 
Ga = (l/v^, ±l/%/2) is a unit vector orthogonal to the 
face of the MBZ where the holon is located. Kinetic en- 
ergy of the holon, eQ(p), is quadratic in the momentum 
p and generally speaking it can be anisotropic. However, 
in LSCO the anisotropy is small and we use the isotropic 
approximation ([2]). 

A very important point is that the argument of e in 
Eq. ([4| is a "long" (covariant) momentum.^®. 



V = -iV + ii7 • [n X Vn] . 



(5) 



An even more important point is that the time derivatives 
that stay in the kinetic energy of the fermionic field are 
also "long" (covariant)'^S,j 



= 9t + -(T • [n X n\ 



(6) 



term" — , — if/'l'^V'a • x n], that is crucially important 
for excitation spectrum and hence for stability of the sys- 
tem with respect to quantum fluctuations. 

Generally speaking, there are also quartic in fcrmion 
operators terms in the effective Lagrangian. However, 
these terms are not important at low doping and there- 
fore we disregard them in ^ . 

The effective Lagrangian ^ is valid regardless if the 
n-field is static or dynamic. In other words it does not 
matter if the ground state expectation value of the stag- 
gered field is nonzero, (n) ^ 0, or zero, (n) = 0. The 
only condition for validity of (|4]) is that all dynamic fluc- 
tuations of the n-field are slow, l/r ^ J , where r is the 
typical time-scale of the fluctuations. We will demon- 
strate below that the dimensionless parameter 



A 



2g^ 



plays an important role in the theory. If A < 1, the 
ground state corresponding to the Lagrangian (0]) is the 
coUinear Neel state and it stays collinear at any small 
doping. If 1 < A < 2, the Neel state is unstable at 
arbitrary small doping and the ground state is static or 
dynamic spin spiral. Whether the spin spiral is static 
or dynamic depends on doping. If A > 2, the system is 
unstable with respect to phase separation and hence the 
effective long-wave-length Lagrangian (|4|) is meaningless. 
Thus, 

A < 1 , Neel state 
1 < A < 2 , Spiral state , static or dynamic 

A > 2, Phase separation . (8) 

For LSCO the value is A « 1.3 - 1.5. 

We would like to stress once more that spin and charge 
to some extent are separated in the effective low-energy 
Lagrangian this is why we use the term "holon" in- 
stead of "hole" . The holon carries pseudospin, it carries 
charge, but it does not carry spin in the usual sense. 
However, it is not the full spin-charge separation like in 
ID models. To illustrate this point, it is instructive to 
look at the holon interaction with uniform external mag- 
netic fielda^!^. 
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{B ■ n)ip'l{a ■ n)ipa 



(9) 



Since we only want to stress the spin dynamics this in- 
teraction does not include terms that originate from the 
long derivative with respect to magnetic vector poten- 
tial — iV — iV — I A, describing the interaction of the 
magnetic field with the electric charge. Clearly the inter- 
action ^ is quite unusual and this is what we call "the 
partial spin-charge separation" . The holon does not in- 
teract directly with the staggered magnetic field (neutron 
scattering). 



The covariant time derivatives result in the "Berry phase III. 



CRITERION OF STABILITY OF THE NEEL 
PHASE UNDER DOPING 



One can consider the coupling constant g in the La- 
grangian Q as a formal parameter. It is clear that the 
Neel order must be stable at a sufficiently small g, 



no = (0,0, 1) 



(10) 



In this case the two hole pockets are populated by holons 
with pseudospin "up" and "down" , and hence the Fermi 
momentum (radius of the pocket) is 



Pf 



(11) 



where x is doping. The Lagrangian (j4|) can be split in the 
diagonal and offdiagonal part with respect to transverse 



(7) spin waves n± = n± = (n^ ± iny)/y/2. 
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q,a) 

^"-^ p+q,E+CO 




P,e 

FIG. 1: Magnon-holon vertex, magnon is shown by the dashed 
hne 



r X± . 2 Ps ( ^ ^ 



(12) 



+ 51(2 -V'ie(p)V'aj , 

(t] - - {p, no[c?n_L X <?]} + \/2.gno[(ea • Vn^) x ct] j 



Here {..., ...} stands for the anticommutator. Using the second quantizatfon representatfon for the n-field, 



with the magnon creation and annihilation operators rnlj_ q and TO±,q, we find the "bare" magnon dispersion 



and the pseudospin flip magnon-holon vertex shown in Fig. [H 

"2 



X-L 



M ^ i 



|\/25(e, . q) + I + 1 [e(p) - e(p + q)]| 



(13) 



(14) 



Looking at (|13p. one can conclude superficially that magnons are hardened by doping. However, they are not hardened, 
they are softened. To see this we need to calculate the magnon polarization operator that is due to Ci. The operator 
reads 

V(,,n\ 2^fn f {V2g(e„-q) + | + l[e(p)-£(p + q)]}' 
/-^Ar(w,q) = Vjp(l-Jp+q) 1—^, 1 ; N I -n 

X± ^ e(p) + w - e(p + q) + 

^J_Y-fn_f , {V2g(e„ ■ q) + I + ^[^(p - q) - e(p)]}' 
X± ^ e(p) - w - e(p - q) + lO 

where /p is the usual Fermi-Dirac step function. Eq. (jlSp can be transformed to 



(15) 



7'Ar(w,q) 

'Po(w,q) 



-g^ + 27'o(c^,q) , 



Then the magnon Green's function reads 



G 



e(p) + ti; - e(p + q) + zO e(p) - cj - e(p + q) + iO 



Xl^ 



-u? - r{uj, q) +iO _ ^2^2 _ 2Vo{uj, q) + ?0 



(16) 



(17) 



I 

The condition of stability of the ground state is the w-axis. Hence this condition is c^q^ > — 27^0(0, q) 
absence of poles of the Green's function at imaginary 
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Xc'^q^, at q <^ pf- Doping x does not appear in this 
criterion. Thus, as it is stated in ([5]), the Neel state is 
stable at any smah doping if A < 1. This criterion has 
been discussed many times, see e.gM>. We have rederived 
it here just to demonstrate how the effective action tech- 
nique works in the known situation. 



IV. THE SPIRAL GROUND STATE IN THE 
MEAN-FIELD APPROXIMATION 

At A > 1 the minimum energy is reahzed with the 
coplanar spiral 



no = (cos Q • r, sin Q • r, 0) , 



(18) 



where Q cx (1,0); (0, 1) is directed along the CuO bond. 
To be specific we assume that Q oc (1,0). Due to the 
holon interaction with the spiral the holon band is split 
in two with cr, = ±1, 



2gQ . 



(19) 



In the ground state only the band with Cj, = +1 is popu- 
lated. Therefore, the Fermi momentum, that is the radius 
of the Fermi circle in each pocket, is 



this section we consider in-plane magnetic excitations. 
An in-plane excitation is described by a small deviation 
(f — (p{t,r) from the uniform spiral ground state (fTS]) . 



(cos(Q • r -1- (/?), sin(Q • r + 93), 0) 



(22) 



In the ground state all the holons are in the pseudospin 
state = 1. The in-plane magnons do not change 
pseudospin, therefore in this section we set everywhere 
(T^ = 1. Substituting expression (j22p in the Lagrangian 
(|4|) we once more find the diagonal and offdiagonal parts 
of the Lagrangian 



Co 



(23) 



2 

E 



(3x 



A 



1 /3 

Ci = ^V't^a (V2g(e„- V)¥J--95---{Z,9(p} 



Here 1 = p+Q/2 is shifted momentum and {..., ...} stands 
for anticommutator. Thus, the "bare" magnon disper- 
sion is given by the same Eq. as for the Neel state, 
but the magnon-holon vertex is smaller than (|14p by the 
factor -\/2, 



(20) 

The point p = corresponds to k = (7r/2,±7r/2) in the 
Brillouin zone. According to (|19p the center of the filled 
holon pocket (ctz = 1) is shifted from this point by — |Q, 
and the center of the empty pocket (ct^ = —1) is shifted 
by iQ. Calculation of energy and its minimization with 
respect to Q gives the following value 
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(21) 



M = i 



V2g{ea ■ q) 



2 



(24) 



A calculation similar to that performed in section IIIII for 
the Neel state gives the following Green's function for the 
field ip that describes the in-plane magnon 



Gi' 



X± 



Vo+iO 



(25) 



The ground state energy of the spiral state is below that 
of the Neel state only if A > 1. 



V. THE IN-PLANE MAGNONS IN THE SPIRAL 
STATE 

To analyze the stability of the spiral state one needs 
to go beyond the mean-field approximation and study 
excitations and quantum fluctuations in the system. In 



where Vo is given by Eq. (jl6[) with the Fermi momen- 
tum ()20|1 . At zero frequency and at small q, q ^ pp, 
the polarization operator is equal to 7^0(0, q) = —^c^q^. 
Therefore, the ground state is getting unstable (poles of 
the Green's function at imaginary w-axis) at A > 2. This 
is the instability with respect to phase separation^Si^ and 
it is fatal for the effective long-wave-length Lagrangian 
dH). Thus, the spiral state is stable at 1 < A < 2, see Eq. 
([8]). We also present here an explicit expression for the 
polarization operator 
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Re 'Po{uj,q) 
Im Vo{'^,q) 
Ro 



e{R', ~ Rt) JRI -Rl- JRI - Rl d{Rt - R^) 



1 2 la 
l3qpF , Rl = ^f3q - uj , i?2 = -^Pq 



(26) 



The Fermi momentum pf is given by (|20p . and 0(x) is 
the usual step function. 
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r 



The g-integrated in-plane magnon spectral density 

cPq 



liniuj) = / /i„(w,q) 



(2n) 



(28) 



is plotted in Fig.[2j3 for doping x — 0.025, x — 0.05, and 
X — 0.1. For zero energy the value of the g-integrated 
spectral density is independent of doping and equals to 
Un{0) = 1/(1 -A/2). 

1 

CO 



0.5 




nl , \ , \ , \ , 

0.5 1 1.5 CO 2 

FIG. 2: A: Plots of 2LuIin{uj,q) versus energy for different 
values of momentum q (offsets). Iin{uj,q) is the in- plane 
magnon spectral density (|27|l . The plots are presented for 
doping X = 0.1, and /3 = 2.7, g ~ 1- Values of q axe given in 
units of the incommensurate vector Q, see (|21|l . 
B: g-integrated in-plane spectral density (|28p for dopings 
X = 0.025, 0.05, and 0.1. The parameters of the effective 
Lagrangian are (3 — 2.7, (7 = 1. 



FIG. 3: The magnon dispersion along Q. The parameters 
are x — 0.1, l3 — 2.7, g = 1. The out-of-plane excitation 
for |g| < Q is shown by the solid line and the out-of-plane 
excitation for |g| > Q is shown by the dashed line. The in- 
plane excitation is shown by the dotted line. The quasiparticle 
residue decays very quickly outside of the dome shown by the 
solid line. The quasiparticle residue at point 1 at the top of 
the dome is Z = 0.8 while the quasiparticle residue at point 4 
that is outside of the dome at the same height is just Z = 0.13. 
The quasiparticle residue of the in-plane magnon at the same 
frequency as the dome height (points 2 and 3) is Z = 0.15. 



It is convenient to define the magnon spectral density 

as 

/„(w, q) = -4pJmG'i„(w,g) . (27) 

Plots of 2ujIin{uJ, q) versus ui are presented in Fig. [2j\ for 
different values of momentimr q (offsets). The doping is 
X = 0.1, and /3 — 2.7, g — I. The narrow peak is the 6- 
function broadened "by hands" to fit in the picture size. 
The corresponding quasiparticle residue is rather small, 
say for q = O.IQ in Fig.[21A_ the residue is Z = 0.39 and it 
very quickly dies out at larger values of q. The magnon 
"dissolves" in the particle-hole continuum. 



To calculate the in-plane spectral density that can be 
observed in neutron scattering one needs to shift mo- 
menta. The Hamiltonian describing the interaction of 
the neutron spin with the n-field reads 

i?^ oc S^^ • n = 5f + i + 5%^+) . (29) 

After the substitution of the in-plane excitation ((22)) . the 
above Hamiltonian reads 
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where k is the momentum transfer and Q the momen- 
tum shift due to the spiral ground state. The scattering 
probabihty for unpolarized neutrons is given by 

I„(Lj,k) = i[/™(a>,k-Q) + /,„(cc;,k + Q)] . (30) 

In Fig. [3] we show by dotted Hues the brunches of hnear 
dispersion that correspond to the quasiparticle peak in 
the spectral function /i„(w,q) plotted in Fig. The 
dispersion is very steep, steeper than the bare magnon 
velocity c, and the corresponding intensities are very low. 



us to resolve the problem because the technique accounts 
exactly all the symmetries. For the out-of-plane excita- 



q,a) 



p+Q+q. e+o) 




q, CO 

p+Q-q, e-o) 



t 



FIG. 4: Magon-holon vertexes with pseudospin flip, magnon 
is shown by the dashed line 



VI. THE OUT-OF-PLANE MAGNONS IN THE 
SPIRAL STATE 

Dynamics of out-of-plane magnons are the most com- 
plicated ones. Stability of the spiral state was ques- 
tioned because of the "marginal" character of the out- 
of-plane excitations if considered in semiclassical 1/S- 
approximation^^'^^. The effective action technique allows 



tion let us write the 7?-field as 

n = — n1 cosQ • r, ^1 — sin Q • r,?iz) , 

and substitute this expression into the effective La- 
grangian ([4]). Neglecting cubic and higher order terms 
in n^, we get the diagonal and the offdiagonal parts of 
the Lagrangian 



(31) 



(3x 
4ps 



A 

— cr, H 

2^2 



P + 2^^^ 



g[Qnz - iV2{ea ■ V)nz] + -h^ 



^ {p, e-^'^-'^lQn, - idn,]] ) ^„ + h.c. 



Here (7+ = (7x + i<^y and the bracket {...,...} stands for anticommutator. According to (|3ip . the "bare" magnon 
dispersion in this case is 



The interaction Ci generates the following two pseudospin-flip vertexes shown in Fig. |4] , 

Ma = [(2p + Q + q)-(Q-q) 

Mb = zy^|5[Q + %/2(e„-q)] + |-^[(2p + Q-q).(Q + q) 
Hence, the magnon polarization operator determined by the vertexes reads 



(32) 



(33) 
(34) 



.9(^-911)-! -1(21 + q)-(Q-q) 



.9(Q + 9||) + f -|(21-q)-(Q + q) 



V 



e(l) - e(l + q) + - A + iO 



e(l) - e(l - q) - - A + iO 



Here gn and q±^ are components of momentum parallel and perpendicular to Q, respectively; /i is the Fermi-Dirac 
step function and l = p + Q/2 is the shifted momentum. Eq. ([51)1 can be transformed to 



P(w,q 
2c2 



9% 



(35) 



f3 



Q • (1 + q/2) 



+ 9M 



e(l) -LO-e{l + (i)-A + iO e(l) + u> - e{l + q) - A + iO ^ 
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This form is explicitly symmetric with respect to uj 
magnon Green's function 



Gout — 



= X I 



q V 1 



where 



ReF+ 



A 



Rl 



1 - 



-UJ and q — q. Integration in ([35)) leads to the following 

(36) 



2 2 / 1 
-1 (1-^ 



■{F+ + F^) + iQ 



2 2 



ImF+ = ^lRl-R\^ 



1 - 



A 



2„2 



4g2 6(76 



6g' 



1 - 



2^2 



i?A = A - + , 



Q 



Q 



I RaQ I 'I i \ 

^4 II 



(37) 



Here d{x) is the step function and i?o is defined in ([26l) . 
The function _F'_ is obtained from by the replacement 
Lo — !■ —UJ in _Ra- 

We define the out-of-plane magnon spectral density as 



^«t(w,q) = -4pjm Gout{(^,q) ■ 



(38) 



Plots of 2a;Joiif(w,q) versus uj are presented in Fig. [5]A. 
for different values of momentum (offsets) and q± = 
0. The doping is a; = 0.1. The narrow peak is the 6- 
function broadened "by hands" , the effective width is 
the same as that for in-plane magnons in Fig. [21 The 
corresponding quasiparticle dispersion fiq is plotted in 
Fig.[3]for direction along Q, (i.e. q = q^, q± = 0) and for 
X — 0.1. The part for \q\ < Q is shown by the solid line 
and the part for Ig] > g is shown by the dashed line. We 
do it to stress that the quasiparticle residue decays very 
quickly outside of the dome. Plot of the residue is shown 
in Fig. [S1A_. To illustrate intensities we compare points 
1-4 in Fig. [31 which correspond to different brunches of 
dispersion with the same frequency. The quasiparticle 
residue at point 1 at the top of the dome is Z = 0.8 while 
the quasiparticle residue at point 4 that is outside of the 
dome at the same height is Z = 0.13. The quasiparticle 
residue of the in-plane magnon at the same frequency as 
the dome height (points 2 and 3 in Fig. [3]) is Z = 0.15. 

An analysis of Eq. (|36p gives the following approximate 
formulas for the dispersion of the out-of-plane magnon 
and for the corresponding quasiparticle residue 



/3xg2c2 
4ps 



(39) 



600 
^500 

T 

^400 
3 

^300 
200 
100 



qil=1.6Q 



q,,= l.2Q 



q|l=0.8Q 



q,i=0.4Q 



qii=o 







0.25 



0.5 



100 

I 

^ 75 
50 
25 
0, 








0.25 



0.5 



FIG. 5: A: Plots of 2ujIout{uj,q) versus energy for different 
values of momentum q (offsets). Iout{uj,q) is the out-of-plane 
magnon spectral density (|38|l . The plots are presented for 
doping X — 0.1 and /3 — 2.7, g = 1- Values of q are given in 
units of the incommensurate vector Q, see (|2ip . 
B: g- integrated in-plane spectral density (|4Up for dopings x = 
0.025, 0.05, and 0.1. The parameters are f3 = 2.7, 5 = 1. 



These formulas have very limited region of validity since, 
as we already pointed out, at larger q the magnon dis- 
solves in the particle-hole continuum. At x = 0.1 the 
equation ([35)1 for fiq agrees reasonably well with the re- 
sult of numerical calculation shown in Fig. [31 At the same 



9 




tral densities 



FIG. 6: A: The quasiparticle residue versus momentum for 
the out-of-plane magnon for the direction along the spiral, 
<? = 91 1 1 l-i- = 0. The vertical line shows the momentum 
where the dispersion vanishes. The doping is a; = 0.1. 
B: The static component of n-field versus doping. The pa- 
rameters are <; = 1, /3 = 2.7. 



time the formula (|39p for the quasiparticle residue only 
poorly agrees with numerics shown in Fig.[6lA_. Certainly 
at very small doping Eq. ([5^ is accurate. 

The g-integrated out-of-plane magnon spectral density 



I out (to) = / /out(w,q) 



(2^ 



(40) 



is plotted in Fig. [5|3 for doping x = 0.025, x — 0.05, and 
X = 0.1. It is peaked at energy Ecross corresponding to 
the top of the dome in Fig. [31 Interestingly, the spec- 
tral density decays almost abruptly to its high frequency 
asymptotic value I{uj) ^ 1 as soon as the magnon is dis- 
solved in the particle-hole continuum. The decay of the 
in-plane g-integrated spectral density shown in Fig. [2] is 
not that steep. 



VII. QUANTUM FLUCTUATIONS AND 
QUANTUM PHASE TRANSITION TO THE 
DYNAMIC SPIRAL PHASE (DIRECTIONAL 
NEMATIC) 

Due to in-plane and out-of plane quantum fluctuations 
the static component of the staggered field n is reduced. 



(41) 



Expectation values ((/9^) and (n^) can be expressed in 
terms of Green's function or in terms of q-integrated spec- 



E 

q 



El UW 



du; 1 



duj 
2^ 



/«„(^),(42) 



These expressions must be renormalized by subtraction 
of the ultraviolet-divergent contribution that corresponds 
to the undoped cr-model. The physical meaning of re- 
lations (|¥T|) and (132) is very simple: the reduction of 
static response is transferred to the dynamic response. 
The most important contribution to quantum fluctua- 
tions comes from out-of-plane excitations with momenta 
q Q (X X. To find this contribution we use the Green's 



function Gout ~ 
([55)) . This gives 



^^p- , where Z and f2 are given by Eq. 



< nl >- 



B 



n/2 



dtdip 



(43) 



1 



t) [(1 -f)2 +2tCOs2, 



Thus, the leading term in the quantum fluctuation scales 
as cx y/x. The subleading contribution to the quantum 
fluctuation scales is x. To find it we have performed nu- 
merical integration in Eq. (I42[) using q-integrated spec- 
tral densities lout and calculated in sections [Vj and 
[VH see Fig. [2^ and Fig.[I^. This gives 



1 - 



2^2^p?/2yr3T7A 



B^/x + 2.6x. (44) 



Certainly the coefficient in the subleading x-term de- 
pends on parameters (a rather weak dependence). The 
value 2.6 in ([44l) corresponds to g = 1 and f3 = 2.7. The 
plot of (n) versus doping x at these values of parameters 
is presented in Fig. [BJ3. 

According to Fig. [BJ3 the static component of n van- 
ishes at X = Xc ~ 0.11. This is a quantum critical point 
for transition to the dynamic spiral. In this phase there 
isn't a spontaneous direction of the n-field, {ft) — 0, but 
the spiral direction (1,0) or (0,1) is still spontaneously 
selected. In our opinion, this is the "nematic phase" 
observed in Refi^. Clearly the value Xc ~ 0.11 is an 
approximate value. In doing the spin- wave theory we as- 
sume that {nl) <C 1, but then, to find the critical 
point we extend this consideration to (n^) ^ 1. This ex- 
tension brings some uncertainty in the value of Xc- We 
also would like to note that the value of Xc is rather sen- 
sitive to parameters. The main sensitivity comes from 
•\/l — 1/A in the denominator in Eq. ([44]) . The value of 
A given by Eq. ([7|) is closely related to the value of the in- 
commensurate vector Q given by Eq. (|2ip . Our estimate 
of Xc is valid for LSCO. 

At X > Xc the spin-wave pseudogap is opened. To 
describe the gapped phase we use the Takaliashi ap- 
proach*^-, see also RefJ^. Idea of this approach is to 
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impose constraint (n) = using the Lagrange multiplier 
method. So we introduce an additional term in the effec- 
tive Lagrangian 



5C = x^K 



1 9 



1 . 

2 - 



(45) 



where is technically the Lagrange multiplier, and 
physically this is the spin-wave pseudogap. The value 
of A, must be determined from the condition 



1 



1 







(46) 



The in-plane quantum fluctuation ((^^) is only very 
weakly (quadratically) dependent on the pseudogap A^. 
The out-of-plane fluctuation (n^) contains a term that 
depends on As linearly. The term comes from the ^Jx- 
contribution in Eqs. To account for the pseu- 




0.5 



0.75 



FIG. 7: The g-integrated magnon spectral density / = -f 
Jout in the gapped "nematic" phase for x = 0.13 and Xc = 
0.11. 

dogap one needs to replace the expression in square 
brackets under the square root in B, see Eq. P5)) . by 



[{l-tf + 2tcos^ (j)] -f 



45 



c2/3(l - IIX)Q' 



-A 



(47) 



A simple calculation with parameters g — \ and (3 = 
2.7 shows that the condition results in the following 
pseudOgap 



As « 2.5(2: — Xc) ■ 



(48) 



This formula is valid only very close to the critical point. 
In this problem one cannot expect a high accuracy from 
the Takahashi-like approach. Therefore, the slope 2.5 
in Eq. (j48p is rather approximate. Finally, in Fig. [7| we 
present the plot of the q-integrated magnon spectral den- 
sity /i„(w) -I- /oMt(w) for X = 0.13 and Xc = 0.11. The 
figure clearly demonstrates that A^ is a pseudogap since 
there is some spectral weight at w < A^. 



given by equation ([2T|) . It depends on the coupling con- 
stant g. Fit of experimental incommensurability! gives 
g « 1 and this agrees remarkably well with prediction of 
the t-t' -t" - J model. 

An important dynamical parameter is E^ross which is 
the height of the dome in Fig. [31 This parameter has been 
systematically studied very recently in inelastic neutron 
scattering^. The experimental values are presented in 
Table I. As soon as the coupling constant g is found from 



X 

EcrossijaeV) 


0.025 0.04 0.05 

7 tt 15 tl 20 tl 


0.07 0.1 

23 40 tl 


phase 
spiral direction 
theory 


insulator 
diagonal 
Refs.3^'3« 


superconductor 
parallel 
present work 



TABLE I: LSCO: Experimental values^ of E^ross versus dop- 
ing X. 

the experimental incommensurability Q we can fit E^^Qsg. 
As we already pointed out above, the present theory is 
applicable to LSCO at a; > Xper ~ 0.055. According to 
Eq. 



l(3xC 



4:PS 



(49) 



Comparing this formula with data at x = 0.07 and x = 
0.1 in Table I we find that /? = 2.65(1 ± 0.1). This value 
agrees reasonably well with the value /3 w 2.2 that follows 
from the t — t' — t" — J model. Note that in principle the 
inverse mass /3 can be somewhat dependent on doping. 
However, the data with error bars are quite consistent 
with x-independent /3. 

Let us discuss also the data at 0.02 < x < Xper = 0.055 
that is relevant to the insulating phase with diagonal 
disordered spin spiral. The corresponding theory has 
been developed in Refs."^^''^*. The incommensurability 
in this case is Q = \/2gx/ps. To fit the experimental 
incommensurability we need g ~ 0.7. This is somewhat 
smaller than the value in the conducting phase. We be- 
lieve that the reduction of g is due to interaction with 
phonons. The point is that g — Zt, where Z is the quasi- 
hole residue. Interaction with phonons in the insulating 
phase can easily reduce the residue by 20-30%. Stability 
of the disordered spiral in the insulating phase is due to 
localization of holes. The Ecross in this case is^^ 



Er 



(50) 



VIII. DISCUSSION AND COMPARISON WITH 
EXPERIMENT 

There are several points that can be directly compared 
with experiment. The incommensurability vector Q is 



where k is the inverse localization length. It is worth 
noting that Eq. 150]) has been derived in Ref.— assuming 
that the binding energy of a hole trapped by Sr ion is 
larger than the magnon energy. The binding energy is 
about 10-15meV. Therefore, strictly speaking, Eq. (I5D1) is 
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applicable only at x = 0.025 since at larger x the energy 
Ecross is getting too big. Nevertheless, we can try to 
apply ([50| to the data at x < 0.055. Fitting the data from 
Table I we find values of k, x = 0.025: k = 0.55 ± 0.2, 
X = 0.04: K = 0.65 ± 0.2, and x = 0.05: k = 0.75 ± 0.2. 
So, there is a hint for a weak doping dependence of the 
inverse localization length k. Most likely the dependence 
is just an imitation of the binding energy correction to 
formula ((50)l . On the other hand, a weak dependence of 
the localization length on doping is quite possible. The 
above values agree reasonably well with the value k « 
0.4 that follows from the analysis of the variable range 
hopping conductivity at a very small doping [x — 0.002), 
see Ref.— ''. 

Near x — 0.12 certain La-based materials in LTT phase 
develop a strongly enhanced static incommensurate mag- 
netic order accompanied by a small lattice deformation 
at the second order harmonic a^i^"'^^ , see also KeiJ^ for a 
review. The measured static magnetic moment ~ O-l/x^ is 
substantially larger than the value that follows from the 
present theory (the unity in the vertical scale in Fig. 
corresponds to the magnetic moment 0.6/ib). There are 
also experimental indications that the spin structure in 
this case is close to collinear— . We strongly believe that 
physics of these materials is somehow related to mech- 
anisms considered in the present paper. On the other 
hand, it is clear that in this case there are some addi- 
tional effects that are not accounted for by the present 
theory. 

The present theory qualitatively explains the direc- 
tional nematic state discovered in underdoped YBCO at 
doping X « 0.08^. For a quantitative comparison one 
needs to analyse the two layer situation. This analysis 
has to include an explanation of a smaller incommen- 
surability compared to that observed in the single layer 
LSCO. 

In the present work we did not account for the su- 



perconducting pairing. The point is that at low doping 
the pairing practically does not influence magnetic exci- 
tations. A different question is how the spiral and the 
corresponding magnetic excitations influence the super- 
conducting pairing. The spin-wave exchange mechanism 
for pairing of holons was suggested in Refs4^i^. The 
mechanism is always working as soon as a short range 
antiferromagnetic order exists in the system. So the su- 
perconductivity peacefully coexists with spin spirals^. 
Moreover, we understand now that the pairing in the 
spiral state is strongly enhanced by closeness to the Neel 
state instability driven by the parameter A. The enhance- 
ment will be considered elsewhere. 

In conclusion, using the low-energy effective fleld the- 
ory we have considered the 2D t-J model in the limit of 
small doping. Quantitatively this consideration is rele- 
vant to underdoped single layer cuprates. We have de- 
rived the incommensurate spin structure (static and/or 
dynamic), calculated spectra of magnetic excitations 
(Figs, m [21 [5]), and considered the quantum phase tran- 
sition to the directional nematic spin-liquid phase. The 
spin wave pseudogap is opened in the spin-liquid phase, 
the q-integrated spectral density in this case is shown in 
FigH 
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